Magnetic ordering, electronic structure and magnetic anisotropy energy in the high-spin Mnio 

single molecule magnet 



O 
O 
(N 

u , 

< 

(N 
(N 



C/3 



I 

O 

o 



> 

o 
o 



I 

O 

o 



X 
S3 



Jens Kortusj^] 

Max-Planck-lnstitut fiir Festkorperforschung, Heisenbergstr. 1, D-70569 Stuttgart, Germany 

Tunna Baruah, N. Bernstein, and Mark R. Pederson 

Center for Computational Materials Science, Naval Research Laboratory, Washington, DC 20375-5000, USA 

(Dated: February 1, 2008) 

We report the electronic structure and magnetic ordering of the single molecule magnet [Mnio04(2,2'- 
biphenoxide)4Bri2]*~ based on first-principles all-electron density-functional calculations. We find that two 
of the ten core Mn atoms are coupled antiferromagnetically to the remaining eight, resulting in a ferrimagnetic 
ground state with total spin S = 13. The calculated magnetic anisotropy barrier is found to be 9 K in good 
agreement with experiment. The presence of the Br anions impact the electronic structure and therefore the mag- 
netic properties of the 10 Mn atoms. However, the electric field due to the negative charges has no significant 
effect on the magnetic anisotropy. 



The interest in magnetic molecular clusters of transition 
metal ions has been continuously growing since the obser- 
vation of magnetic bi-stajbility of a purely molecular origin 
in the so-called Mni2-ac,l!l which shows a magnetic hystere- 
sis cycle below 4 K similar to that observed for bulk mag- 
netic materials. The magnetic bi-stability associated with the 
hysteresis cycle has created an interest in these clusters for 
information storage, although at low temperature quantum ef- 
fects affect the rejszersal of the magnetization, resulting in steps 
in the hysteresis.ElThis phenomenon of quantum tunneling of 
magnetization is governed by the magnetic anisotropy energy 
(MAE)H barrier which is due to directional dependencies of 
the spin-orbit-coupling operator. 

Recently, Pederson and Khanna have developed a raethod 
for accounting for second-order anisotropy energies. □ This 
method relies on an exact and simple expression for the spin- 
orbit coupling operator used in a second-order perturbative 
treatment to determine the dependence of the total energy 
on spin projection. Initial applications to the uniaxial Mni2- 
ac lead-to a density-functional-based second-order anisotropy 
energytl of 55.J K, in agreement with the experimentally de- 
duced valuesBD of 54.8(3) K or 55.6 K. 

Because the second-order anisotropy energy scales with the 
square of the magnetization it was generally believed that 
a high-spin ground state S would be beneficial for a large 
barrier. The [Mnio04(2,2'-biphenoxide)4Bri2]^^ cluster has. 
been reported to have n S — 12 high-spin ground stateO 
but only a small energy barrier of about 7.7 K. In this work 
we investigate the electronic and magnetic properties and the 
magnetic anisotropy energy of this high-spin single molecule 
magnet. The information obtained here may be useful in the 
search for single molecule magnets with a greater magnetic 
anisotropy. 

Fig. |l] shows the structure of the [Mnio04(2,2'- 
biphenoxide)4Bri2]*^ molecular magnet. The 10 Mn atoms 
form a tetrahedron like structure with Mn atoms at the cor- 
ners and at the middle of the tetrahedron edges. Two of 
the Mn atoms, the top and the bottom spheres in Fig. |l], are 
coupled antiferromagnetically to the rest of the Mn atoms. 
The Mn atoms are bridged by O atoms. The magnetic 




FIG. 1 : The tetrahedron like structure of the ten Mn atoms and the 
surrounding organic rings. The Br atoms are not displayed for clarity. 
Density isosurfaces for 0.03 e/a% for majority (dark) and minority 
(light) spins on Mn atoms are shown. The plot clearly shows that 
the magnetic moment is localized at the Mn atoms, and it directly 
confirms the antiferromagnetic coupling of two of Mn atoms (large 
light spheres) to the remaining Mn atoms (large dark spheres). 



core is further stabilized by organic rings that are also con- 
nected to the O atoms. The negatively charged cluster is 
compensated by [(CH3CH2)3NH]2[Mn(CH3CN)4(H20)2] in 
the molecular crystal, but experimental results suggest that 
the magnetic anisotropy is due to the localized valences of 
the 10 Mn atoms.Q In order to use the higher symmetry 
and make the problem computationally feasible, we carried 
out our calculation on the negatively charged [Mnio04(2,2'- 
biphenoxide)4Bri2]^~ cluster which contains 1 14 atoms. The 
eight symmetry operations reduce the complete cluster to 1 8 
inequivalent atoms. 

The theoretical studies were caiTied out using a linear com- 
bination of atomic-orbitals molecular-orbital (LCA£UMO) 
approach within the density-functional frameworkMlJ The 
molecular orbitals were expanded as linear combinations of 
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TABLE I: The Gaussian basis set used for the calculation. The min- 
imum and maximum exponent a of the bare Gaussians, the number 
of bare Gaussians, the number of contracted s-, p- and d-like basis 
functions. 
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Gaussian functions centered at the atomic sites. The calcu- 
lations were carried out at the all-electron level and the mul- 
ticenter integrals required in the solution of the Kohn-Sham 
equation were .calculated by integrating numerically over a 
mesh of points.liil 

Our density functional-based calculations were performed 
with the all-electron Gaussian-orbital based Naval Re- 
search L|abejcatorAU|Molecular Orbital Library (NRLMOL) 
program,l!ill3l!iLjOli3 using the Perdew-Burke-Ernzerhof 
(PBE) generalized-gradieuL approximation for the exchange 
and correlation functionalO NRLMOL combines large Gaus- 
sian orbital basis sets, numerically precise variational integra- 
tion and an analytic solution of Poisson's equation to accu- 
rately determine the self-consistent potentials, secular matrix, 
total energies and Hellmann-Feynman-Pulay forces. The ex- 
ponents for the single Gaussian have been fully optimized for 
DFT calculations. The basis set for the Mnio cluster con- 
sisted of a total of 3756 contracted orbitals. The minimum and 
maximum exponent of the bare Gaussians, the number of bare 
Gaussians, and the number of contracted s-, p- and rf-like ba- 
sis functions are given in Table | for each atomic species. The 
contraction coefficients for atomic orbitals were obtained by 
performing an SCF-LDA calculation on the spherical unpolar- 
ized atom where the total energy of the atom was converged to 
within 10 meV. The basis functions that do not correspond to 
atomic wavefunction were constructed from the longest range 
bare Gaussians in the basis set. 

Here we repeat some of the formulas needed for discussion 
of the magnetic anisotropy energy. The same definitions and 
notation are used as in Ref.^. In the absence of a magnetic 
field the second-order MAE A2 resulting from the spin-orbit 
coupling, for an arbitrary symmetry, reduces to 

aa' ij 

which is a generalization of Eq. (19) of Ref.^ The matrix el- 
ements = {x'^\Si\x'^ ) implicitly depend on two angles 
(6, (3) defining the axis of quantization. The matrix elements 
Mfj^ , which are related to of the induced orbital moment, are 
given by 

j^jaa' ^ {(l>la\Vi\4>ka'){(l)ka'\Vj\(l)la) 
kl ^''^ ~ 

where (j)i„, (p^a and ei^ are, respectively, the occupied, unoc- 
cupied and the corresponding energies of states. Vi is same 



as defined in Eq. (7) of Ref.^ and is related to derivatives of 
the Coulomb potential. The matrix elements can be evaluated 
by integrating products of the Coulomb potential with par- 
tial derivatives of the basis functions. This procedure avoids 
the time consuming task of calculating the gradient of the 
Coulomb potential directly. 

In addition to the magnetically interesting complex, the 
crystal also contains single Mn complexes to balance the. 
charges. Using high-field EPR spectroscopy Barra et aln 
found that this [(CH3CH2)3NH]2[Mn(CH3CN)4(H20)2] unit 
is paramagnetic. We also find that this unit is paramagnetic, 
with the Mn atom in a +2 charge state and a spin of S=5/2. 
The complex exhibits easy-plane behavior with an energy well 
of 0. 1 K. We therefore focus the remainder of our work on the 
Mnio unit only. 

Calculations on a 5*= 12 high spin state revealed that this 
spin state would not be magnetically stable because the Fermi 
levels in the majority and minority spin channel would not be 
aligned. The Fermi level misalignment indicated further trans- 
port of electrons of the minority to the majority spin channel. 
As a result we obtained a 5*= 13 high spin state as the mag- 
netic ground state instead of the 5"= 12 state obtained from 
high-field EPR spectroscopy. Our result is consistent with 
experiment since it is difficulllo differentiate experimentally 
between the two possibilities.t3 

A plot of the spin density shown in Fig. |l] clearly confirms 
the antiferromagnetic coupling of two of the Mn atoms to the 
remaining Mn atoms. The spin density around the minority 
spin Mn and the four majority spin Mn at the corners of the 
tetrahedron-like magnetic core show a nearly spherical spin 
density which one would expect for a closed (f' shell. The 
other four Mn atoms which are on the edges of the tetrahe- 
dron (large dark spheres) show a less spherical spin density, 
indicating another charge state for these atoms. In order to 
analyze the magnetic ordering we calculated the spin density 
in spheres around the atoms. This gives a measure of the local- 
ized moment at the atom, but will generally underestimate the 
exact value. Using spheres with a radius of 2.23 as around 
the Mn atoms we obtained for the three non-equivalent Mn 
atoms (majority spin tetrahedron edge, majority spin tetrahe- 
dron vertex, and minority spin tetrahedron edge) a local mo- 
ment of 3.61,4.33 and -4.33 /is, respectively. This result sug- 
gests an ionic picture that the first Mn has an Mn'^+ {S = 2) 
state, whereas the other two are Mn^+ (S — 5/2). This pic- 
ture is fully in accord with the spin density plot in Fig. |l]. Due 
to the symmetry of the cluster, the two types of majority spin 
Mn atoms have a multiplicity of 4 whereas the minority spin 
Mn atom has a multiplicity of 2, resulting in the previously 
mentioned S'=:4x2 + 4x5/2 --2x5/2 = 13 magnetic 
ground state. 

The electronic density of states (DOS) for the majority and 
minority spin channels is shown in Fig. |[ For each spin, the 
DOS is further decomposed into the 3d contributions of all 
Mn atoms. It is evident from the plot that the states around 
the Fermi level are clearly connected with 3d states of the Mn 
atoms. This result also agrees well with the experimental pic- 
ture that the states near the Fermi level are well localized and 
do not show strong hybridization with other atoms, although 
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FIG. 2: Electronic density of states (DOS) broadened by 0.54 eV of 
Mnio in the spin S = 13 configuration. For each spin the total DOS 
and the projected DOS of all Mn(3d) are presented. The vertical line 
indicates the Fermi level. 



we find some oxygen, bromine and nitrogen contributions for 
the occupied states. 

Starting from the experimental geometryllJ we carried out 
about 30 steps of a conjugate-gradient algorithm using the 
Hellmann-Feynman-Pulay forces for optimization of the ge- 
ometry. For each new geometry we calculated the complete 
Hamiltonian of the magnetic anisotropy and the second-order 
contribution to magnetic anisotropy barrier DS^. In accord 
with experimental data, we find that the Mnio single molecule 
magnet is an easy-axis system. The baiTier showed no strong 
dependence on geometry varying between 8.8 K and 10.4 K 
with a value of 9.5 K for the lowest energy geometry. Ex- 
pressing the barrier in the form of a simple spin Hamiltonian 
H — DSl we obtained a value of D = —0.056 K. For the cal- 
culation of the spin orbit matrix elements, we included all va- 
lence electrons and all unoccupied states in an energy window 
of 13.6 eV above the highest occupied state. The difference 
between the second-order treatment and exact diagonalization 
of the Hamiltonian including the spin-orbit matrix elements, 
which includes some higher order effects too, was of the order 
of 0.1 K. 

Eq. ^ shows that the baiTier is related to matrix elements 
between occupied and unoccupied orbitals in the majority and 
minority spin channels. In order to give a deeper insight into 
which states are forming the barrier, we analyze these contri- 
butions in more detail. First, we focus on the contributions of 
the different spin channels. Table [| summarizes the result in 
form of the D parameter allowing only a given spin channel, 
for example including only matrix elements Mf-^ between 
occupied majority states and unoccupied minority states, in 
the calculation of the barrier. All matrix elements from the oc- 
cupied majority electrons prefer an easy-axis system, whereas 
the matrix elements from the occupied minority spin channel 
would result in an easy-plane system. Only due to the larger 
values of the contributions of the occupied majority spin chan- 
nel the system ends up as an easy-axis system. This destruc- 
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TABLE II: The contributions of the different spin channels (see 
Eq. ^ to the magnetic anisotropy parameter D and the magnetic 
anisotropy energy DS^. 



tive interference between the different spin channels seems 
to be the reason for the small barrier compared with Mni2. 
In Mni2 constructive interference between the different spin 
channels was observed. 

Besides the spin channel contribution, we can analyze 
which electronic states contribute most to the matrix elements 
K'^-^ . In Fig. |3]we display plots of the square of the wave- 
functions of the occupied majority state and the unoccupied 
minority state which contribute to the matrix element Mff 
with the largest absolute value. The view is from the top with 
respect to Fig. [l|. It is clearly visible that the states of inter- 
est are d-states localized at the same Mn atom. In this case, 
the states are localized at the minority spin Mn atoms (light 
spheres in Fig. p. In order to emphasize the d character of 
the wavefunctions, we have chosen the top view, although the 
wavefunctions of the other minority Mn atom are just below 
the top ones and are not visible. 

While Mn is the only magnetically active species in the 
complex, the remaining atoms affect the magnetic properties 
of the molecule. In particular the electric field of the twelve 
Br~ ions can affect the MAE through its effects on the elec- 
tronic structure and on the spin-orbit coupling. Effects of the 
Br ions on the electronic structure are of a chemical nature, 
and a detailed analysis is beyond the scope of this work. How- 
ever, direct effects on the spin-orbit coupling energy could 
raise the possibility that even small variations in the posi- 
tions of these ions could affect the magnetic properties of the 
molecule. 

To measure the effects of the Br^ ions on the electronic 
structure we redid the calculations with various subsets of the 
Br~ ions. For each removed Br we also removed an extra 
electron, keeping the remaining molecule isoelectronic with 
the original complex. For these systems with either zero, four, 
or eight Br atoms, we observed a range of behaviors. In some 
cases the electronic structure near the Fermi level was similar 
to the original molecule, although it never showed truly rigid- 
band-like behavior In other cases, however, the electronic 
structure changed significantly, sometimes completely clos- 
ing the HOMO-LUMO gap. Associated with these electronic 
structure changes were large changes in the MAE, including 
changes in the magnitude of the anisotropy barrier, as well as 
instances of changes to an easy-plane system. 

To measure the direct effect of the electric field of the Br 
ions on the spin-orbit coupling, external Coulomb potentials 
which acted to cancel the long range affects of the Br an- 
ions were added and the spin-orbit interaction and magnetic 
anisotropy were recomputed. This neutralized the electric 
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FIG. 3: Isolines at 0.005 e/a| of the square of the wavefunctions (occupied majority state (a) and unoccupied minority state (b) ) which 
contribute most to the matrix elements M^f . The view is from top with respect to the earlier figures. It is clearly visible that the matrix 
element connects majority and minority d-states at the same Mn atom. 



field due to the Br^ ions near the Mn sites without changing 
the electronic structure of the molecule. We tested the effects 
with neutralizing charge distributions of various widths, and 
by neutralizing four, eight, or all twelve Br anions. The MAE 
changed by less than 1 K in all of these calculations. We there- 
fore conclude that the electric fields created by the Br^ ions 
do not have a significant effect on the magnetic properties of 
the molecule. 

In conclusion, we present a study of the electronic and mag- 
netic properties of the Mnio single molecule magnet. We 
confirm the experimentally suggested magnetic ordering, al- 
though we find that a state with S — 13 is the magnetic ground 
state in contrast to the S = \2 state suggested from high field 
EPR measurementsJa In agreement with experiment we find 
the Mnio unit is an easy-axis system with a small barrier of 
9.5 K and the compensating cluster in the molecular crystal. 



which has one Mn atom with S* — 5/2, is an easy plane sys- 
tem with a MAE of 0. 1 K, negligible compared with the Mnio 
unit. We show that the magnetic anisotropy is determined 
by a competition between different spin channels involved. 
The electric field caused by the negative charges of Br anions 
has no significant direct effect on the spin-orbit coupling or 
the MAE, although their chemical interactions do have sig- 
nificant effect on the electronic structure and therefore on the 
MAE. The states most important for the magnetic anisotropy 
energy involve transitions between occupied majority and un- 
occupied minority d-states at the same Mn-atom. 
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